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^ '. Abstract 

The effects of the gravitational back reaction of cosmological perturbations are investigated in 
a cosmological model where the universe is dominated by phantom energy. We assume a COBE 
^vq ■ normalized spectrum of cosmological fluctuations at the present time and calculate the effective 

I/-} . energy-momentum tensor of the gravitational back-reactions of cosmological perturbations whose 

> ; 

wavelengths at the time when the back-reactions are evaluated are larger than the Hubble radius. 
' Our results reveal that the effects of gravitational back-reactions will counteract that of phantom 

© : energy sooner or later and can become large enough to terminate the phantom dominated phase 

before the big rip as the universe evolves. This arises because the phase space of infrared modes 
Oh' grows very rapidly as we come close to the big rip. 
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I. INTRODUCTION 



The analysis of data from supernovaeQ], CMbQ and WMApJjJ strongly indicates the 
existence of dark energy which dominates the present universe and_drives the acce_lerat 
ing cosmic expansion. Many models, such as cosmo" 
essence 0], braneworldj?]], Chaplygin gas[^, quintom 
to explain the dark energy. A phantom field 
tion of state(u> = p/p < —1) in contrast to quintessence energy {w > —1) or cosmology 
constant (w = — 1), appears as another possible candidate of dark energy and has received 



■ iical constant quintessence Q], k- 
and holography[lO|, are proposed 
which has the super- negative equa- 



increased attention recently. It has an unusual kinetic term in its Lagrangian which gives rise 
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to some strange properties, such as the violation of the dominant energy condition 
and the increase of its energy density with time. As bizarre as it may appear, such terms 
could arise in a variety of theories 14 1 . 

It can be shown that once our universe enters the phantom energy dominated phase, the 
scalar factor will blow up in a finite proper time due to excessive expansion. This arises 



because 
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|15| the energy density of phantom fields increases with time instead of red- 
shifting away as the matter or radiation energy densities or as the energy density of ordinary 
quintessence. Thus the increasing phantom energy will ultimately strip apart gravitationally 
bound bodies and cause a cosmic doomsday or big rip 15]. The above conclusions are based 
upon a constant negative value of w. However, if the value of w could change during the 
evolution of the universe and then in principle the big rip can be avoided [12]. Attempts have 
also been made toward avoiding the big rip by modifying the original Caldwell's phantom 
model HQ, 0,0,0- 

Note that it has been argued that a singularity can also develop 
at a finite future time even if p + 3p is positive 0]. 

Here we would like to point out that the big rip may be avoided in phantom cosmology 
even with a constant negative value of w without adding new physics. We will demon- 
strate that the gravitational back reaction of cosmological perturbations may terminate the 
phantom dominance before the big rip occurs. 
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II. GRAVITATIONAL BACK REACTION 



Considering the fluctuations of metric and matter and expanding the Einstein equation 
to second order, due to non-linearity of Einstein equation, the second order equation can not 
be satisfied and the back reaction of these fluctuations to the background must be considered 
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which is characterized by a gauge-invariant effective energy-momentum tensor r^ iv 

r = (T (2) —G {2) ) (1) 

(2) (2) 

where T^J and G)J express the second order metric and matter perturbations and pointed 
brackets stand for spatial averaging. This formalism can be applied to both scalar and tensor 
perturbations and applies independent of the wavelength of the perturbations. The effects 
of gravitational back reactions have been studied in the context of cosmological models, 
where, for example, they have been used to address the issue of dynamical relaxation of the 
cosmological constant [3| and the possible termination of quintessence phase Here we 
will discuss the back reactions in phantom cosmology. 

We will assume that w is constant and w < — 1, and the universe is not phantom domi- 
nated until the time t m . That is, t m is the time of equal phantom energy density and matter 
energy densities. The universe is matter dominated if t < t m and is phantom dominated 
if otherwise. In the following we discuss the back reaction of cosmological perturbations 
during the phantom dominated phase (t > t m ). However, if phantom exists all the time, the 
back reaction should too before t m . Here we assume that the back reaction is negligible or 
its effects are reinforcing before the phantom dominated era , since if otherwise the phantom 
phase could not occur. As the amplitude of each fluctuation mode is small, we need a very 
large phase space of modes in order to produce any interesting effects. During the phantom 
dominated phase, both the scalar factor a and expansion rate H increase with time, so Hub- 
ble distance decreases and the phase space of infrared modes grows. Hence we expect that 
the effects of the back reaction of infrared modes will grow and we therefore only focus on 
the back reaction of infrared modes on the evolution of the universe dominated by phantom 
energy. 

A simple effective action of a phantom field can be expressed as 

c phan = - l -{d^) 2 -v{v), (2) 
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where V(tp) is the potential of the phantom field. Thus we can obtain the energy momentum 
tensor 



^d a <pd a ip + V{ip) 



(3) 



In longitudinal gauge the metric with scalar perturbations can be written as 



ds 2 = a 2 (r])[(l + 2$)V - (1 - 2$)&ijdx x d3p] 



(4) 



where $ is the Bardeen potential. In addition to the geometrical perturbations, one must 
also consider the perturbations of the phantom field 5<p during the phantom dominated era. 
Expanding the Einstein equation to the first order in $ and Sip, we obtain the gauge-invariant 
equations of motion for small perturbations 



$' + h$ = -A-n:G(p' 5(p 



(5) 



(6) 



$" + 3h& + {2ti + h 2 )<S> = -4ttG[-<p%$ + (p' V + V, v a 2 5ip] , 



(7) 



where a' represents the derivative with respect to conformal time r\ and h = a' /a. Eqs. fl5l6l7|) 
come from the Oi, 00 and ij components of the Einstein equation respectively. Subtracting 
Eq. (jSJ) from Eq. (JZJ) and using Eq. (jSJ) and the equation of motion for the phantom field, 
we get a second order partial differential equation for $ 



$" _ V 2 $ + 2 ( h - ^ $' + 2 [ti - ] $ = . 

<PoJ V <Po, 



(8) 



This equation is the same as that in which the matter perturbations are induced by a usual 
scalar field 



Thus, for long-wavelength perturbations we have 



H . , 
1 / adt 

a 



Here Ak is an integration constant. For short- wavelength perturbations one gets 



(9) 



(10) 



where a dot denotes the derivative with respect to coordinate time t. 
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Expansion of the energy-momentum tensor T^ v and Einstein tensor G^ v to the second 
order in <3> and Sip yields the non vanishing components of the effective back reaction energy- 
momentum tensor 



too = g^[12# <<M>> " 3($ 2 ) + 9a- 2 (V$ 2 )] - \{5^) 
-ia- 2 ((VM 2 ) + \v, w (V) + 2V, V (Sty) , 



(11) 



and 



Tij ^'StG 



(24# 2 + 16H)($ 2 ) + 24#($$) + ($ 2 ) + 4($$) - -a~ 2 (V$ 2 ) 



-i\/, w (V) + 2^($^)} , (12) 
where H = a/a. 

Let us now apply these equations to the phantom dominated era. For t > t m the solution 
for the scale factor is given by 



a(t) = a(t m ) 



w + (1 + w) — 

tm. 



3(1+to) 



(13) 



Apparently a diverges when t = wt m /(l + w) = ti rip . During the phantom dominated era 
the energy density evolves as p p han ~ a~ 3( - 1+w \ thus p p han{t) is related to p(t m ) by 



Pphani^) 



[-W + (1 + w)t/t r 



(14) 



Here p(t m ) ~ l/(67rGr^). The energy density p p han{t) increases with time and diverges when 
t = t brip . Using p = wp, p = — |y? 2 — V((p) and p = — |y? 2 + U(y?), we obtain 



and 



Vltp) 



■ 2 = ~(W + l)p(t m ) 

V " [-w + (l + w)t/t„ 



(-w + l)p(t m ) (-w + l)p{t m ) 



(15) 



2[-w+ (l + w)t/t m Y 



exp 



-2(l + u;)c 



t m y/-(w + l)p(t m ) 



(16) 
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Combining Eq. © and Eq. ([130. we get in the long- wavelength limit 

f 0A k = A k , 

*fc = i (17) 
I A fe [l - lnx(t)] = A fc z(£), iu = -5/3, 

where /3 is an integration constant and x(t) = | — As the universe evolves from t m to 
tbrip i x{t) varies from 1 to and z(t) from 1 to oo. It then follows that = $^ = for 
w 7^ —5/3 and 

®k = \ ; : $fc , (18) 

3t mX (t)z(t) *' 1 J 

and 

•■4 

for u> = —5/3 . Defining p br = Tq and pb r = and using Eqs. ()5I15I16I18I19[) . the 

expression for t^ v can be simplified to 

^ F <$ 2 >, W ^ -5/3 , 



'67rGt2J-uH-(l+™)t/*mp 



(20) 



-Pfer, w 7^ -5/3 , 

/v: ' 2 ^ 2 2W w = -5/3. m] 



%*GV mX (tY ^ *(t) 

Here the subscript 6r stands for back reaction. So the equation of state parameter for the 
dominant infrared contribution to the back reaction is given by 

-/ ^-5/3, 

^ ~ 1 2+3,(t)-2 2 (Q' r /o ^ 

k l+3,(t)-2,(t) 2 ' ^ ~ °/ t3 • 

Hence, for iu 7^ —5/3, the contribution of infrared modes to the energy momentum tensor 
which describes the back- reaction takes the form of a negative cosmological constant, whose 
absolute value chane; function of time. The result of this case is similar to that 



23, 



24j . which has been used to address 



obtained in an inflationary background cosmology 

the issue of dynamical relaxation of the cosmological constant [25]. Meanwhile, for ir 
—5/3, Wb r = —3/2 at t — t m , and w^ r grows as time goes on and approaches to —1 at tb r i p . 
Using Eq. ffHll2"Uj). we have 

-(l-w)($ 2 ), w ^-5/3 



Pbr 

P^n"l|(^ + 4-2)($ 2 ), «. = -V3 



(23) 



If the above ratio is negative, the phantom energy will be counteracted by the effects of the 
back reaction. When the ratio becomes negative unit, the phantom phase will end and the 
universe will re-enter the matter dominated era. In order to determine the value of this ratio, 
it is pivotal to evaluate the two-point function ($ 2 ), which can be obtained by integrating 
over all Fourier modes of $: 

where hi = ciiHi and kt = a(t)H(t) are infrared and ultraviolet cutoffs respectively. The 
infrared cutoff can be chosen as the length scale above which there are no significant fluctu- 
ations. If we assume that there was a period of inflation in the early history of our universe, 
we could take the infrared cutoff as the Hubble radius at the beginning of the inflation. The 
ultraviolet cutoff will be taken as the Hubble radius at time t when the strength of the back 
reaction is to be evaluated. 



III. TERMINATION OF PHANTOM DOMINATED PHASE 

We assume that our universe now is dominated by the phantom energy, and for simplicity, 
t m ~ ^0; which means that the phantom just begins to dominate at today. As a matter of 
fact, one of the main features that distinguish the phantom energy from a cosmo logical 
constant or quintessence is that the onset of phantom energy dominance happens at the 
very last moment. So our assumption is fairly reasonable. Now we want to relate <3> p fc 
defined as the value of <3>fc during the phantom dominated phase, and $ m fc, which is the 
corresponding value just before the phantom energy begins to dominate. Let us consider 
separately modes whose wavelengths are outside the Hubble radius now and modes whose 
wavelengths are inside the Hubble radius today and exit the Hubble radius before time 
t(t > t ). For the 
Bardeen potential 



brmer case, <& p k(t) and $ m fc(io) can be related by the conservation of the 
2llQ 



C=-~ - + 25 

S 3 1+w V ; 

In the matter dominated phase, since w = and $ m fc = 0, we obtain C(^mJfc) = 5$ m fc/3. For 
w 7^ —5/3, during the phantom dominated era, = 0. So we have ((§ p k) = (5 + 3w)/(3 + 
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3w)$ p fc. It follows from the conservation of ( that 

$ = 5(1 + w) $ mfc . (26) 
p 3w + 5 K J 

For modes whose their wavelengths are inside the Hubble radius today and exit the Hubble 
radius at time tn given by a(tn)H = k. From Eqs. (jl0ll5ll7|) . we know how $ evolves 
with time on the wavelength scale smaller and larger than the Hubble radius. Thus we can 
express $ p fc(t > t H ) in terms of t , t H and <& m k{to) 

$ , (f) = <$> mk (t )- ^ r . (27) 

pkK ' mfeV 0J [-wt + (l + w)t H ] 

For the case of w = —5/3, we have in the matter dominated era, from Eq. (jHJ), that $ m fc = 
in the long-wavelength limit. For modes whose wavelengths are outside the Hubble 
radius now, using Eq. ()17|). we have 

%k{t) = -$mk(t )z(t) . (28) 

For those modes whose wavelengths are inside the Hubble radius today and exit the Hubble 
radius before time t(t > t ), using Eqs. (jl0ll5ll7|) . we obtain 

**w = *^Jsi§)' (29) 

We now assume that the spectrum of cosmological fluctuations at the present time is 
normalized by the recent observations of CMB anisotropies and the fluctuational spectrum 
has the form 

P(k) = ^|^(to)| 2 = C (t^—) X , (30) 
2vr 2 \kcoBEj 

where kcoBE = <ya H and C 1 ^ 2 ~ 10 -5 . By the joint analysis of the Maxima-1, Boomerang 
and COBE cosmic microwave anisotropy results and the estimate of the systematic errors, 
the blue spectrum tilt should be in the region of < A < 0.27 j^J 1 . In this case a ~ 7.5. 
Substituting Eqs. (I2h127l3()|) into Eq. flU), we obtain, w ^ -5/3 



2 25C(w + If f ko A _ x Ct 2 Q f kt k x - x 

{<P[t} } ~ (3w + 5f(aa H o y J ki + J^H& J ko [-wt + (1 + w)t H ]* dk ' 



1 According to the more recent WMAP data the upper bound could be slightly larger [EI ■ 



The first integration in the above can be evaluated as 

25C(w + l) 2 rko 

9i ~ 



.A-l 



(3w + 5) 2 (aa H ) x j, 
25C(w + l) 2 
(3w + 5) 2 h ' 



dk 



(32) 



where we have defined 



fi 



ln^fl. 



(33) 



A = , 

A^a-^l - ^(a F )- A ], A + . 

For an e-folding number as large as 60 ~ 70, hx(aoHo/ki) ~ O(10 l ). Therefore, since 
C ~ 10~ 10 , gi <C 1. The second integration is given by 

k^ 1 



92 



cti 



(aa H v ^ x 



or Jk [-wt + (1 + w)t H ] 
(l + 3w) C 



-dk 



6(H-iu) + (l + 3w)Aa 



w + (1 + w) 



t 



l+3w 



A-2 



(34) 



Apparently g% ~ when t ~ to and it increases with time and approaches infinity at t — tb r i P - 
Thus hereafter, we will discard g\ in discussing if and when the effects of back reaction can 
terminate the phantom dominated phase before everything is torn apart at the big rip. 
For w = -5/3, the substitution of Eqs. ()2SI29l3()jl into Eq. flMJ) leads to 

25Cz(t) 2 f k ° \ i Czitf f kt k^ 1 
1 k x ^dk + 



^COBE Jk 

25 



k COBEJk z(t H ) 2 x(tH) 2 



dk 



2Cz{tf 



fi + f*(t) 



(35) 



where 



^COBE 



„A-1 



fco z(t H ) 2 x(t H y 



-dk 



2a-\l + A)e- 2(1+A) ^[2(l + X)z(t)} 



x(t) 



-2(1+A) 



(36) 



a x z{t) 

Here Ei(x) is the exponential integral function. Appealing to the fact that when x — > oo 
Ei(x) can be expanded as 



Eiix 



~ e x (x~ 1 + x 2 + 



we have that 



/ 2 (t) ~ a- A X (t) 



-2(1+A) 



Z (ty 



2(1 + A) 



+ ■ ■ 



(37) 



(38) 
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when t approaches Urip- The above result blows up at the big rip, so we will discard in 
Eq. flUSJ). Using Eqs. J23EU33E5) , we obtain 



_ (l+3w)(l-t») C 
Pbr J 6(l+w)+(l+3w)\ a x 



-w + {l + w)± 



l+3-w 



A-2 



"3(1+™)" 1 , , _ ,„ 

- 1 , w -5/3 , 



/V """ 1 ^[l + 3^(t)-2z(t) 2 ]/ 2 (t), w = -5/3. 



(39) 



3 

It is interesting to note that when w ^ —5/3, the ratio is always negative as long as w < — 1 
and becomes negative unity before t — > h rip . While, for the case of w — —5/3, it can be 
shown that the ratio is positive at present and remains so for a period of time. The plot 
of this ratio vs t/to in Fig. ^ shows that it turns negative approximately at t ~ 1.81£q. 
This indicates that the back reaction reinforces the phantom energy in an early period of 
phantom dominated universe and then counteracts it. At the big rip this ratio becomes 
negative infinity. Therefore, the phantom dominated phase will be terminated sooner or 
later before the big rip by the back reaction effects. The behaviors of the energy density 
of the back reaction and that of the phantom background as a function of t/to are plotted 
in Fig. 121 and Fig. |H1 There we can see that the ratio of \p p han/ Pbr\ reaches unity (i.e., 
Logi (\p p han/ Pbr\) becomes zero) before the big rip and thus the phantom phase terminates. 
Note that the termination point appears very close to the big rip in the figures. However, 
since to(~15Gyr) a is very large number, the actual time interval measured in terms of years 
is not small (see also the Table). Let us now discuss in more detail when this happens. For 
the case of w ^ —5/3 the time when the phantom phase is terminated by the effects of 
gravitational back reaction, i.e., when Pbr ~ —1, can be explicitly given as 

Pphan 



t 



3(l+m) 



/ [6(1 + w) + (1 + 3w)Xa x 1 (i+3»)A+6(i +u ,) \ 



1 +w 

Define t' as the time of termination of the phantom phase before the big rip, we have 

t 



t — tbrip t 



1+W 



3(1+™) 

6(l + it;) + (1 + 3w)Xa x t ^ 



(1+3iu)A+6(1+uj) 

(41) 



(l + 3u;)(l-tw) C 

For a fixed w, t' is an increasing function of A. This is similar to the case of quintessence 
in that the effect of back reaction is proportion to the blue tilt. For scale invariant spectra, 
i.e., A = , the effect of back reaction for quintessence is negligible In contrast, for the 
phantom case, it is easy to see from Eq. ()41j) that the effect of back reaction can still become 
large enough to terminate the phantom phase. When w = —5/3, it is almost impossible to 
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TABLE I: The time of termination of the phantom phase before the big rip: to = 15Gyr, a = 
7.5, C = 10- 10 . 



t'(yr)(X = 0) t'(yr)(X = 0.27) 



-1.1 4.2 x 10 b 6.4 x 10 8 

-1.3 9.6 x 10 5 1.9 x 10 7 

-1.5 5.1 x 10 5 5.2 x 10 6 

-5/3 3.4 x 10 5 2.7 x 10 6 

-1.8 3 x 10 5 2 x 10 6 

-2 2.4 x 10 5 1.3 x 10 6 







Q- 
Q. 



-2 

-4 • • 

-6 



. . w=-5/3 
— w=-1.5 



1 1.2 1.4 1.6 1 . 

t/t 



FIG. 1: Pbr/Pphan vs i/io is plotted for ui = 1.5 and w = —5/3 at the early epoch of phantom 
dominated phase with A = 0.27, a = 7.5, C = 10~ 10 . The vertical coordinate has been scaled by 
C = 10" 10 . 

solve Eq.(jnHl) to get an analytical expression for the time when Pbr/Pphan — ~ 1> but we can 
resort to numerical techniques. Listed in the Table are results of the calculations for the 
time of termination of the phantom dominated phase for both w ^ —5/3 and w = —5/3. 



IV. DISCUSSIONS 



The main conclusion we can draw from the Table is that the big rip can be avoided by 
the gravitational back reaction of cosmological perturbations. A comparison of our results 
with that of Ref. 3| for the case of w — — 1.5 shows that our Solar system could be saved!. 
We can also see from the table that for a fixed w the greater the blue tilt, A, the sooner the 
termination of the phantom phase. However, for a fixed A, the greater the w, the earlier the 
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Q. 1- 





Big Rip: 
\ 

\ 



60 

o 

J -1 w=-1.5 



Termination 



2.99 3 

t/t 



FIG. 2: Plotted is Log w {\ Pphan / Phr \) vs t/t for w = -1.5 with A = 0.27, a = 7.5, C = 1(T 10 . The 
step length of the horizontal axis is 0.0005. Note that to ~15Gyr 
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Q_ 



» 

o 

J 1 









B~icj Rip^ 














w=-5/3 


Termination 



2.495 2.5 

t/t 



FIG. 3: Log w (\p p han/pbr\) vs t/t is plotted for w = -5/3 with A = 0.27, a = 7.5, C = 10~ 10 . The 
step length of the horizontal axis is 0.0002. Note that to ~15Gyr 

termination of phantom phase before the big rip. The physical reason is that the greater 
the w, the farther away we are from the big rip and thus there is more time for the infrared 
modes to accumulate. Finally, let us consider a case of pure theoretical significance, i.e., the 
case in which A — > oo, then we find t' ~ 13Gyr for w ^ —5/3 and t = 15Gyr . This means 
that it is not possible to kill the phantom dominated phase the time just when it kicks in. 

In summary, We have, assuming a COBE normalized spectrum of cosmological fluctua- 
tions at the present time, calculated the gravitational back-reaction effects of cosmological 
perturbations whose wavelengths at the time when the back-reactions are evaluated are 
larger than the Hubble radius. Our results reveal that the gravitational back-reactions are 
growing with time and could become large enough to terminate the phantom dominated 
phase before the big rip occurs. An interesting feature to be noted of the gravitational 
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back reactions is that their effective energy momentum tensor is that of some form of "mat- 
ter" which has ne gati ve energy density and positive pressure. This form of "matter" was 
postulated in Ref. jl6[ 
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